Scale-dependent bias of galaxies and fi-type distortion of the cosmic microwave 
background spectrum from single-field inflation with a modified initial state 
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We investigate the phenomenological consequences of a modification of the initial state of quantum 
fluctuations of a single inflationary field. While single-field inflation with the standard Bunch-Davies 
initial vacuum state does not generally produce a measurable three-point function (bispectrum) in 
the so-called squeezed triangle configuration (where one wavenumber, k, is much smaller than the 
other two, k <C ki ~ £2), allowing for a non-standard initial state produces an exception. Here, we 
calculate the signature of an initial state modification in single-field slow-roll inflation as it would 
appear in both the scale-dependent bias of the large-scale structure (LSS) and /x-type distortion 
in the black-body spectrum of the cosmic microwave background (CMB). We parametrize the ini- 
tial state modifications and identify certain choices of parameters as natural, though we also note 
some fine-tuned choices that can yield a larger bispectrum. In both cases, we observe a distinc- 
tive k~ 3 signature in LSS (as opposed to the k~ 2 of the so-called local-form bispectrum). As a 
non-zero bispectrum in the squeezed configuration correlates one long-wavelength mode with two 
short- wavelength modes, it induces a correlation between the CMB temperature anisotropy observed 
on large scales with the temperature-anisotropy-squared on very small scales; this correlation per- 
sists as the small-scale anisotropy-squared is processed into the ^i-type distortion of the black-body 
spectrum. While the correlation induced by the local-form bispectrum turns out to be too small to 
detect in near future, a modified initial vacuum state enhances the signal by a large factor owing to 
an extra factor of ki/k compared to the local form. For example, a proposed absolutely-calibrated 
experiment, PIXIE, is expected to detect this correlation with a signal-to-noise ratio greater than 10, 
for an occupation number of about 0.5 in the observable modes. Relatively-calibrated experiments 
such as Planck and LiteBIRD should also be able to measure this effect, provided that the relative 
calibration between different frequencies meets the required precision. Our study suggests that the 
CMB anisotropy, the distortion of the CMB black-body spectrum, and the large-scale structure of 
the universe offer new ways to probe the initial state of quantum fluctuations. 



I. INTRODUCTION 

While cosmologists have accumulated extensive evi- 
dence for an early-universe inflationary period, the cause 
and dynamical specifics of that epoch remain unclear. 
Current and upcoming measurements will provide in- 
creasingly precise measurements of the effects of infla- 
tion, demanding that theorists persist in relating these 
observations to inflation's underlying mechanism. Pri- 
mordial non-Gaussianity is a popular discriminant among 
the proposed models of inflation (e.g., [1-3]). 

The scalar curvature perturbation, which appears 
in the space-space part of the metric in a suitable gauge 
as gij = a 2 (t)e 2 ^ 5ij (where a(t) is the Robertson- Walker 
scale factor) , is a convenient quantity relating the observ- 
ables such as the cosmic microwave background (CMB) 
and the large-scale structure (LSS) of the universe to the 
primordial perturbations generated during inflation. In 
particular, this quantity is conserved outside the horizon 
for single-field inflation (e.g., [4]). We shall define the 



two-point function (power spectrum, denoted as P^(fc)) 
and the three-point function (bispectrum, denoted as 
B^(ki, fc 2 , k)) of £ in Fourier space as follows: 

(CkXk) = (2w) 3 5(k 1+ k)P c (k), (1) 
(C kl Ck 2 Ck) =(2vr) 3 ( 5(k 1 +k 2 +k) J B c (fc 1 ,fc 2 ,fc). (2) 

The current data constrain the shape of P((k) as P^(fc) oc 
fc" s " 4 with n s = 0.96 ± 0.01 [5, 6]. 

The so-called local-form bispectrum defined as [7, 8] 

^ ocal (fci,fc 2 ,fc) = ^/ NL [P c (fc!)P c (fc 2 ) + (2 perm.)] , 

is particularly interesting, both because a detection of the 
primordial bispectrum at the level of /nl 3> 1 would dis- 
favor single-field inflation [9-12] and because it is easy to 
measure the primordial signal since few late-time effects 
can produce the local-form bispectrum. The most impor- 
tant contamination of /nl known to date is due to the 
lensing-Integrated Sachs- Wolfe (ISW) effect bispectrum 



2 



[13], which can be calculated precisely and removed. The 
contamination of /nl due to non-linearity in the photon- 
baryon fluid has been shown to be at most one [14]. 

The local-form bispectrum has the largest signal in the 
so-called "squeezed triangle configuration," for which one 
of the wavenumbers, say, k, is much smaller than the 
other two, k <C k\ « k 2 . This can be seen from Eq. (3): 
as P^(k) cx k~ 3 for a scale-invariant spectrum (n s = 1), 
the bispectrum is maximized when k is taken to be small. 
In this limit, one finds: 



B local (fci,fci,fc -> 0) -+ ^-f^P^P^k) cx 



1 



kfk 3 ' 



(4) 



for a scale-invariant spectrum. 

Recently, Agullo and Parker have shown that a non- 
standard initial state of quantum fluctuations generated 
during single-field inflation can enhance the bispectrum 
in the squeezed configuration by a factor of ki/k, i.e., 
-> 0) cx l/(fc?fc 4 ) [15]. This would have pro- 
found implications for observations of the bispectrum in 
the squeezed configuration. For example, the signature 
in the bispectrum of CMB of this model was investigated 
in a paper by one of the authors [16], who found that the 
model could produce a measurable local /nl signal in the 
CMB. 

The primordial bispectrum in the squeezed configura- 
tion was initially constrained mostly by measurements of 
the temperature anisotropy of the CMB [5, 17]. How- 
ever, over time, tools for observing the bispectrum have 
proliferated, providing a variety of ways to compare in- 
flationary models. In this paper, we will explore two such 
methods: 

1. In the large scale-structure (LSS) of the universe, 
the local-form bispectrum leaves a signature by 
contributing a scale-dependence to the halo bias, 
b(k) [18-20]. For the local-form bispectrum, the 
scale dependence goes as 1/fc 2 ; however, for a mod- 
ified initial state, this scale dependence can become 
1/k 3 . 

2. Anisotropy in the so-called /x-type distortions of the 
black-body spectrum of the CMB can be correlated 
with the CMB temperature anisotropy measured 
on large scales. This correlation can be used to 
measure the bispectrum in the squeezed configura- 
tion but with a larger value of k\ jk than previously 
thought possible [21]. 

This paper is organized as follows. In Section II, we 
review the model under consideration. In Section III, 
we give the form of the bispectrum and comment on po- 
tential uncertainties in the results. In Section IV, we 
discuss a useful approximation to the bispectrum in the 
squeezed configuration. In Section V, we calculate the 
signal of this model in the scale-dependent bias of LSS. 
In Section VI, we calculate the signal of this model in the 
/z-type distortion of the CMB black-body spectrum, cor- 
related with the CMB temperature anisotropy on large 
scales. Finally, we conclude in Section VII. 



Throughout this paper, we shall set M v \ = l/y/SirG = 
1, and use the cosmological parameters given by the 
WMAP 5-year best-fit parameters (WMAP+BAO+i7 
ML; [22]): fl M = 0.277, fl A = 0.723, h = 0.702, 
n s = 0.962, and A 2 ( {k = 0.002 Mpc^ 1 ) = 2.46 x 10~ 9 , 
unless stated otherwise. 



II. ACTION AND MODE FUNCTION 

We consider here single-field slow-roll inflation with a 
canonical kinetic term, where the action (to lowest order 
in slow-roll) can be written as [9] 

S =S 2 + S3 , 



HCd- 2 (. 



(5) 



We expand the curvature perturbation into creation, a^, 
and annihilation, ak, operators (not to be confused with 
the Robertson- Walker scale factor, a(t)): 



C(*,t) 



d 3 k 

(2tt) 3 



a k u fe (i)e lk - x + aK(*)e- 1 • (6) 



Usually, one chooses an initial state so that a comov- 
ing observer in the approximately de Sitter spacetime 
observes no particles (i.e., for this observer dk|0) = 0). 
This implies a positive-frequency mode function given by 



u k {r l ) = ^ r ^={l + ikr l )e- ikr ', 



(7) 



where n = J* ^4pr is the conformal time; for future ref- 



erence, we note 



du k _ H 2 
dn 6 



qe 



-ikr] 



(8) 



While this is certainly a reasonable assumption, it is 
an assumption, and all assumptions must be tested by 
observations. Thus, a responsible scientist should ask: 
"If the initial state of ( was not in this preferred vacuum 
state (known as the Bunch-Davies state), what are the 
implications for observations?" Our goal in this paper 
is not to construct candidate models of a modified ini- 
tial state, but to study phenomenological consequences 
of such a modification, i.e., to let our observations tell us 
about the initial state of quantum fluctuations. 

Once we adopt this approach, the next question is: 
"How should we parametrize a modified initial state?" 
We will represent a modified initial state as a Bogoli- 
ubov transformation of the above Bunch-Davies mode 
function: 



Uk(v) = ot k u k (n) + f3 k u* k (r]). 



(9) 
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This is not the most general form one can write down 
(see, e.g., [23]), but it provides us with a reasonable start- 
ing point. In line with our previous goal, we will take the 
Bogoliubov coefficients as given rather than trying to de- 
rive them from a fundamental theory. From the com- 
mutation relation of creation and annihilation operators, 
the coefficients at and (3k must satisfy \ak\ 2 — \flk\ 2 = 1- 
We also find that the occupation number of particles Nk, 
i.e., the expected number density of particles with mo- 
mentum k, is given by \(3k\ 2 - 

These Bogoliubov coefficients, ak and (3k, encode in- 
formation about physics on scales where we have lim- 
ited information; thus, they can vary widely without in- 
consistency However, we can place some constraints on 
the coefficients by demanding that the theory reproduce 
the observed power spectrum (including the spectral tilt, 
n s = 0.96 ± 0.01) and that the energy in the fluctuations 
not back-react on the background inflaton dynamics [23- 
25] . These requirements can be satisfied in a fairly natu- 
ral way if we suppose that the coefficients are such that 
(Nk) « Nk$e~ k / fe <™*, where the cutoff momentum fc cut 
must be specified. The values allowed for Nk,o depend on 
the value of fc cut [25]; for fc cut ~ y/M p \H, i.e., the scale 
of inflation, Nk,o can be of order unity. Additionally, if 
we suppose that the smallest primordial scales observable 
today come from momenta sufficiently smaller than fc cu t> 
then (Nk) ~ N k _ = N, i.e., roughly constant in k. Re- 
membering that (Nk) — |/3fe| 2 and that only the relative 
phase between ak and (3k is significant, we parametrize 



a k 



= s/T+Nj 



Pk^VN. (10) 



There is still uncertainty with respect to 6k- As ex- 
plained further in [16], we identify two scenarios as plau- 
sible behaviors: 1) 9 k ~ krj , where fc|r7 | ~> 1 for relevant 
k, and 2) 9k ~ const = 9. In the latter scenario, one can 
tune the value of 9 to give larger effects; we will gener- 
ally show results that assume the value of 9 that gives 
the largest signal. In this sense (and for another reason 
discussed in [16]), we consider the former scenario to be 
more conservative. 



III. POWER SPECTRUM AND BISPECTRUM 

The power spectrum of £ on super-horizon scales, kr\ <C 
1, which seeds the observed fluctuations, is given simply 
by P c (k) = \ti k (v^0)\ 2 [16], i.e., 



*(*) = fr^l«* + &l 



which becomes (using Eq. (10)) 
H 4 1 



(11) 



Pc(k) = ^2p( 1 + 2N + 2VN(N + 1] C0S9k ) ■ (12) 
The calculation of the bispectrum requires more 



thought. Formally, it is given by [9] 

(C kl WCk 2 WCk(t)) 

= -*/ dt'([( kl (t)C k2 (t)C k (t),H int (t')}), (13) 

J to 

where the interaction Hamiltonian, H mi , is given by 
/ dt'H int (t') = -S 3 with #3 given by Eq. (5). We would 
then specify the initial state at the initial time, to, or 
equivalently at the initial conformal time, rj a . For the 
action given by Eq. (5), one finds 



B c (k u k 2 ,k) = 2i 



1 2 ' 2 ' 1 — 



k' 1 h 2 h' 2 
1 2 



-U kl U k2 Uk 



'/() 



U h U h U h +C.C. 



7\3 "fci "fe 2 a k 



(14) 



In this paper, dots will denote derivatives with respect 
to t and primes will denote derivatives with respect to r\. 
For the standard calculation, we take the Bunch-Davies 
initial vacuum state, given by ak — 1 and (3k = 0, for 
all modes into the infinite past, rjo — > — oo (i.e., t — > 0). 
For this case, there is an accepted prescription for calcu- 
lations: we take r]o — > r/o + ie|r/o| ; giving r\ an imaginary 
component when its absolute value is large [9] . The expo- 
nential terms in the integrand like e l ( kl+k2+k * >Vo (see (8) 
for their origin) would ordinarily oscillate rapidly at very 
early times but are suppressed by the imaginary part of 
jjo- Note that this suppression depends on ki + k 2 + k > 0. 

However, when we allow for a more general initial state, 
we can have (3k ^ resulting in terms like e 4 (-' £ i+ fe 2+' c )»jo ; 
e i(-k 1 -k 2 +k)r] ; e t c Furthermore, one may object (e.g., 
for reasons of renormalizability) to setting initial condi- 
tions in the infinite past, especially if some of the modes 
are excited (i.e., (3k ^ 0); instead, one might prefer that 
initial conditions be set at some finite time. If we ig- 
nore this objection for a moment, one can still suppose 
that r]o — > i]o + ie\i]o\. By triangle inequalities (e.g., 
k\ < k 2 + k, etc.), the exponentials are still suppressed 
except at the precise folded limit k\ = k 2 + k (note that 
this would result in Eq. (15) but without the exponen- 
tials). 

In this paper, however, we will generally take the objec- 
tion seriously and suppose that initial conditions were not 
set infinitely far in the past. Unfortunately, this draws us 
into an area of active research which does not offer a def- 
inite formalism for calculations. Here, as in [26] (though 
see [27, 28]), we will adopt the "Boundary Effective Field 
Theory" approach to non-Bunch Davies initial conditions 
[29] , which like the other available approaches is not with- 
out problems (e.g., [30]). In this approach, one cuts off 
the integral given in Eq. (13) at a finite 770, where the 
initial conditions are set. 

We shall assume that for excited modes k (i.e., where 
(3k ^ 0), k\t] \ » 1 so that k was deep inside the horizon 
at the initial time. This can be explained as expressing 
the requirement that the process for mode excitation was 
causal and thus, could only excite subhorizon modes. 
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Performing the integral given in Eq. (14), one obtains 
the bispectrum [16] 



B ( (ki,k 2 ,k) = 



IH^kjk 2 + k 2 k 2 + k 2 k 2 

2 r/,2 



x 3? 



1 



k\ + k 2 + k 



-F n 



Z->3 Z->3 7^3 
1 2 

^ - e i(fei+fe 2 +fc)»7o^ 



+ k l + k 2 -k ta ^V ) 
+ k,-k 2 + k apa \ J 



-fci + k 2 + k 



- {3aa 



where 



F XYZ = (a kl + (3 kl )(a k2 + (3 k2 )(a k + h)X* k J k \Z* k 

-("fa +Pk 1 )( a t 2 + P*k 2 )( a l + P k ) X kl Y k2 Z k > 

(16) 



for = /3, /3° = a; Fxyz gives information about 
the initial conditions at r/o- Note that we ignore a field 
redefinition term (derived in [9]) that is negligibly small 
for the purposes of this paper. 

First, note that we recover the standard Bunch-Davies 
result [9] if we set F aaa = 1, F aa(3 = F a(3a = F (3aa = , 
and rjo — > — (1 — ie)oo. In the squeezed limit, k <C fci m k 2 



and we get B c 2^P c (k 1 )P c (k), where ^ w 1CT 2 
is the slow-roll parameter, which is equivalent to ,/nl = 
0(1CP 2 ) (sec Eq. (4)). If we restore the field-redefinition 
piece we ignored, we obtain the full standard squeezed- 
limit bispectrum: Bq — > (1 — n s )P^(ki)P^(k) [9]. 

Since we have assumed fc|?7o| 3> 1, the exponentials in 
the bispectrum (15) oscillate rapidly and can, to a decent 
approximation, be ignored. Then, one sees that the bis- 
pectrum peaks in the so-called "folded triangle configu- 
ration," where one of the wavenumbers is approximately 
equal to the sum of the other two, i.e., k ~ fci + k 2 , 
k 2 ~ ki + k, or k\ k 2 + fc; this was noted earlier by 
[26, 27, 31]. Since the local bispectrum has no corre- 
sponding peak, this regime provides a way to distinguish 
the shape of this bispectrum from a purely local form. 
We shall come back to this point in Section VII. 

We can also investigate the squeezed configuration 
k <C ki s=y k 2 ; this configuration is in fact a special case 
of the folded limit fci ~ k 2 + k when we additionally sup- 
pose that k is much smaller than fci or k 2 . In this limit, 
the third and fourth terms are larger than the first and 
second by a factor of fci/fc ^> 1; the bispectrum becomes 

12 > 

B( oc 8^ 2 -- k L P(;(ki)P^(k) (with a proportionality factor 



- e tkr i°)]). Note that 



\a kl + fi kl \- 2 \a k + /3 k \- 2 $t[F af3a (l 

this is enhanced relative to the local form in the squeezed 
configuration [15]. 

We should highight that the exponential terms cannot 
be completely ignored [27] because they prevent the bis- 
pectrum from blowing up in the folded limit. In particu- 
lar, the factor _ k + k +k [1 - e t( -~ kl+k2+k ^' no }, which seems 



to blow up in the folded limit if one ignores the expo- 
nential, actually goes as — ir/ + 0((—ki + k 2 + fc)r/o). 
Accounting for this behavior plays a role in the usefull- 
ness of the approximation we demonstrate in the next 
section. 



IV. APPROXIMATION TO THE BISPECTRUM 
IN THE SQUEEZED CONFIGURATION 

While the full form of the bispectrum given by Eq. (15) 
is complicated, the observables that we shall discuss in 
this paper (the scale-dependent halo bias in LSS and the 
anisotropy in the /u-type distortion of the CMB black- 
body spectrum) depend primarily on the squeezed config- 
uration, k <C fci s=y k 2 . Therefore, it is useful to find an ac- 
curate approximation to the bispectrum in the squeezed 
configuration. 

In [16], the author expanded to the lowest order in k 
(fc here is equal to k$ in [16]) after averaging over the 
exponential. Specifically, he approximated 



^ _ e i(k\-k 2 + k)i]a ^ 



(17) 



This result is also consistent with a prescription of ignor- 
ing oscillating terms by taking 770 — > r/o + «e|?7o| for large 
1 770 1, as discussed in the previous section. 

When we do not ignore oscillating terms, the approx- 
imation demonstrates the correct scaling on large scales 
but it is off by a factor. This arises because the approxi- 
mation does not properly account for the oscillatory be- 
havior of Eq. (17) at small k. 

Fortunately, we can come up with a better approxima- 
tion. Observe that, when calculating observables, the bis- 
pectrum is usually multiplied by a function and then inte- 
grated over some of the wavenumbers (see, e.g., Eq. (32) 
below). Let us focus on the integral over k 2 . Note that 
the limits of integration for k 2 are k 2 € [\ki — k\, k\ + fc]; 
in the squeezed limit, the function multiplying the bis- 
pectrum will vary little over this small range, while the 
oscillatory terms like the left hand side of Eq. (17) will 
vary very rapidly. Thus, we can perform the k 2 integral 
only over the rapidly oscillating term, e.g., 



■ki+k ^ _ gi(fci-fc 2 +fc)?7o 

dki I dk 2 • • ■ 

ki —k 



fci — k 2 + k 



I 



dki 



ki+k ^ _ e i(k 1 -k2+k)rj Q 

k 2 =k! \ J kl ^ k k\—k 2 + k 



(18) 



For this integral, we find 



fel+fe l _ e i(ki-k 2 + k)r]o 

fei-fe ki-k 2 + k 



= [7 - Ci(-2fc??o) + log(-2fc??o)] + i Si(-2fc77o) 
« 7 - Ci(-2fcr? ) + log(-2fc77o) , 



■5 



where 7 « 0.5772 is Euler's constant, Ci(z) = 
— J z °° dtcos(t)/t is the cosine integral, and Si(z) = 
Jq dtsin(t)/t is the sin integral (which is ~ it/2 for 
z > 1); in the last line, we have dropped the second 
term since it becomes increasingly unimportant for large 
km- 

If we perform this new approximation, we find that, for 
9 w kr] , the chief contributor to the squeezed bispectrum 
looks like 



B 



8 k Oikrio 

C,fc«fc 



1(ki,k2,k) 



H 6 



ft k x k 2 k A 
x N(l + N) 



x 1[ 7 - Ci(-2fcr7o) + log(-2fcr/ )]. 

(19) 



For 9k ~ const = 0, we find 
iJ 6 1 



Blf™™\ kl ,k 2 ,k) 



2 /cifc 2 fc 4 
AT(1 + iV)(3-cos26>) 

+ y/N(l + N)(l + 2N) cos 6> 
1 

log(-2fc?7o)]. 

(20) 



x -[ 7 -Ci(-2fc % ) 



These two equations provide useful approximations to the 
bispectrum from a modified initial state in the squeezed 
configuration. 

Note that we can also view this approximation as find- 
ing a sort of average for the left hand side of (17), i.e., 
that 



I _ e i(ki-k 2 +k)r) 



ki - k 2 + k 



1 

2k: 



avg 
ki+k ^ 

dk 2 — 

k\—k 



O i{k 1 -k 2 +k)7} 



k\ - k 2 + k 
_ 7 - Ci(-2fc77 ) + log(-2fcr/ ) 
2k 

Thus, ignoring the oscillating terms in the bispectrum 
(as in [16]) is equivalent to neglecting a factor |[ 7 — 
Ci(— 2kr] ) + log(— 2fc?7o)], so that the regime investigated 
in [16] differs from the one here by a factor of order unity. 



V. SCALE-DEPENDENT BIAS 

How can we measure B^{k\,k 2 ,k) observationally? 
Obvious observables are the bispectrum of the CMB tem- 
perature and polarization anisotropy, and that of the 
matter density distribution in LSS. These observables are 
(in linear theory) related to Bq (k\ , k 2 , k) in a straightfor- 
ward way [8, 32]. 



A much less obvious observable is the power spectrum 
of dark matter halos (in which galaxies and clusters of 
galaxies would be formed) . Dark matter halos are formed 
only at the locations of peaks of the underlying matter 
distribution. While the power spectrum of the underly- 
ing matter distribution is insensitive to the bispectrum, 
the power spectrum of peaks is sensitive to the bispec- 
trum as well as to higher-order correlation functions [33] . 
This leads to a remarkable prediction: one can use the 
observed power spectrum of the distribution of galaxies 
(and of clusters of galaxies) to measure the bispectrum 
of primordial fluctuations [18-20]. 

In general, as the power spectrum of peaks (hence ha- 
los) is different from that of the underlying matter dis- 
tribution, we say that halos are biased tracers of the un- 
derlying matter distribution [34]. The degree of bias is 
often parametrized by the so-called "bias factor," b(k), 
defined as 



b 2 (k) 



-Phalo(fc) 



Prr 



r(*)' 



(21) 



Alternatively, one may define b(k) as the ratio of the 
matter-halo cross power spectrum to the matter power 
spectrum. 

On large scales, where the matter density fluctuations 
are still in the linear regime, b(k) approaches a constant 
for Gaussian matter density fluctuations, b(k) — > b\. 
However, the presence of the primordial bispectrum leads 
to a non-trivial fc-dcpendcncc in b(k), and this is called 
a "scale-dependent bias." 

Building on the previous work on the peak statistics 
[20, 33], Desjacques, Jeong, and Schmidt arrived at the 
following formula for b(k) [35]: 



Ab(k,R) = 2 



M R (k) 



(61 - l)6 c + 



d\nF R {k)/d\nR 
din a r I din R 



(22) 



where Ab(k,R) = b(k,R) - h, S c = 1.686, R is re- 
lated to the mass of halos under consideration as M — 
^-n M p c R 3 , and p c = 2.775 x 10 11 h 2 M Q Mpc~ 3 is the 
present-day critical density of the universe. The various 
functions are defined by 



1 



d 3 fci 



4o%P c (k)J (2^)3 

xB c (fti,|ki + k|,ft) 

4 ee / ^-P d k)M%(k), 



M R (k) 
W R (k) 



J (2tt) 3 
2k 2 D(z) 
5n M H 2 
3ji(kR) 
kR 



Mfl(*i)M«(|ki + k|) 

(23) 
(24) 



T(k)W R (k) , 



(25) 

(26) 
(27) 



where T(k) is the linear transfer function normalized such 
that T(k) — > 1 for k — > 0, and D(z) is the growth factor 



6 



of linear density fluctuations normalized such that (1 + 
z)D(z) — > 1 during the matter era. (For example, D(0) ~ 
0.7646 for Cl M = 0.277, fl A = 0.723, and w = -1.) 

Before we show the numerical calculations, we will first 
try to analytically explore Eq. (22) for the case of a mod- 
ified initial state, allowing us to estimate the bispectrum 
shape. 

An important observation in what follows is that 
W R (k) oscillates rapidly for fc > 1/R, so that the in- 
tegral of F R {k) is dominated by fci ~ 1/R. Therefore, 
when we are interested in fc fci, the integral of .F^(fc) 
is dominated by the squeezed configuration, and is ap- 
proximated as 



1 



4a|P c (fc) 



d 3 h 
(2tt) 3 



A4(fci)B c (fci,fci,fc)- (28) 



We can insert the squeezed configuration local bispec- 
trum and calculate 



F R (k) » -/nl- 
5 



(29) 



The second term in the parenthesis in Eq. (22) vanishes 
for this case and we find Ab(k,R) cx 1/fc 2 [18-20]. 

For the bispectrum for a modified initial state, which 
goes as l/fc 2 fc 4 , we instead find 



h(R) 



(30) 



where 



h(R) = — t" 



d 3 h 
(2tt) 3 



hMlik^Pdh). (31) 



One may interpret fci as a characteristic wavenumber for 
the short-wavelength mode in the squeezed configuration. 
Thus, we expect the modified-state bispectrum to pro- 
duce a scale-dependent bias which grows faster (by a fac- 
tor of k\jk « l/(kR)) for small values of k than that for 
the local-form bispectrum. 

What about the second term in the brackets in 
Eq. (22)? If we note that the extra k factor in the in- 
tegrand of Eq. (31) (as compared with the integrand for 
a 2 R ) is evaluated at roughly 1/R, we get ki(R) ~ 1/R 
and din F R (k)/d\nR ss — 1. 

On the other hand, a R is dominated by the power 
spectrum of matter density fluctuations at k 
1/R. Approximating the power spectrum of mat- 
ter density fluctuations as a power-law near k w 
1/R, i.e., M R {k)P c (k)\ kaR -i cx k n ^ R \ one obtains 
d In a R /d In R = — [n e ff(-R)+3]/2. For example, n e s(R) = 
-2.2, -1.8, and -1.6 for R=l,5, and 10 ft," 1 Mpc (or 
M = 3.2 x 10 11 , 4.0 x 10 13 , and 3.2 x 10 14 ft" 1 M Q ), 
respectively. 

Therefore, while this second term changes the ampli- 
tude of Ab(k, R) by a factor of 1 + [» cff(fl)+ 3 ] 5 c (b 1 -i) » [t 
does not change the fc-dependence of Ab(k, R). We thus 
expect the fc-dependence of the scale-dependent bias for 




0.001 0.01 
* [h Mpc -1 ] 



FIG. 1. Scale-dependent halo-bias from single-field inflation 
with a non-standard initial state, using a smoothing scale of 
R — 1 h~ x Mpc. The occupation number is N = 0.5, the 
slow-roll parameter e = 0.01, and the initial conformal time 
1 770 1 = 1.0 x 10 6 Mpc (the results are insensitive to the exact 
choice of 770, so long as it is large). The bottom three (thicker) 
lines show the more natural case where 8k ~ krjo, while the 
top two (thinner) lines show the case when 8k ~ const is 
chosen to give the maximal halo bias. The dashed lines show 
the new approximations given by Eqs. (19) and (20), while the 
dot-dashed line shows the approximation used in [16] (which 
is equal to Eqs. (19) and (20) without the last lines). 



a modified initial state to be given by Ab(k,R) oc 1/fc 3 . 
This scaling was also predicted by [36] . 

In principle, the second term can change the amplitude 
of Ab(k, R) by a large factor for low-mass halos whose 
bias is closer to unity [35] . Nevertheless, as we are focused 
on the shape of Ab(k, R) rather than on the amplitude, 
we will ignore this factor. Then, Eq. (22) simplifies to 



Ab(k,R) _ 1 5 C 
bi-1 ~8ir 2 <7 2 R M R (k)P c (k)k J Q 



dfci fciA^fl(fci) 



x / dfc 2 k 2 M R {k2)B c (ki,k 2 ,k) 
J\ki-k\ 

1 

207T 2 D{z)a R n M H 2 k 3 P ( (k)T(k)W R (k) 



dfci kfT(k 1 )W R (k 1 







ki+k 



dfc 2 klT(k 2 )W R (k 2 )B (; (k l ,k, 2 ,k), 

(32) 



which agrees with the formula first derived by [20]. Here, 
a R = a R / D(z), which is independent of z. 

To evaluate Eq. (32) we will use R = 1 h^ 1 Mpc 
(corresponding to M = 3.2 x 10 11 h^ 1 M Q ). Also, in 
order to determine the factor H 6 /4> 2 which appears in 
the bispectrum, we use the WMAP 5-year normalization, 
fcgP c (fc )/(27r 2 ) = 2.41 x 10~ 9 for fc = 0.002 Mpc" 1 [22], 
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Ab/(b-l) x D(z) 
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FIG. 2. Scale-dependent bias from the local-form bispectrum 
(dot-dashed line) versus the modified initial-state case de- 
scribed herein (solid line). The parameters here are the same 
as in Fig. 1, with /nl = 1 for the local-form bispectrum. The 
difference in scaling between the models is quite evident. 

in Eq. (12): 

= ^- (i + 2N + 2VmTT)cose) , 

2ir z o7T e V / 

(33) 

where e = (l/2)cf> 2 /H 2 is the slow-roll parameter; for 
8k ~ kr)o, the term in parenthesis simply becomes (1 + 
2N). This relation gives H 2 for a given e, TV, and 9. 

We now insert the full bispectrum (Eq. (15)) into 
Eq. (32) and numerically integrate for the halo bias. 
Fig. 1 shows the results of the numeric integration for 
e = 0.01, N = 0.5, and 9k for both scenarios (9k = kr/a 
and 9 = const, with 9 chosen to maximize Ab). We do 
find the expected k~ 3 scaling, which can also be seen by 
comparison with the local form in Fig. 2. 

Fig. 1 also shows the halo bias as calculated from the 
approximation given in Eq. (19), as well as from the ear- 
lier approximation from [16] (which is Eq. (19) without 
the last line). One sees that the full calculation and the 
new approximation are very similar except on very small 
scales (near the smoothing scale); by contrast, the old 
approximation is off in absolute scale. (There is also a 
slight change in shape, due to the terms that depend on 
kf] ). 

VI. tt-TYPE DISTORTION OF THE 
BLACK-BODY SPECTRUM OF CMB 

A. Motivation and background 

Diffusion damping of acoustic waves heats CMB pho- 
tons and creates spectral distortions in the black-body 
spectrum of the CMB [37]. However, this distortion is 
erased, maintaing a black-body spectrum for the CMB, 



as long as photon non-conserving processes are effec- 
tive. According to [38], double- Compton scattering 
(e~ + 7 — > e~ + 27) is an effective thermalization process 
for z > Zi ~ 2 x 10 6 . After this epoch, however, this pro- 
cess shuts off and the spectral distortions from diffusion 
damping cannot be smoothed from the CMB spectrum. 
Since elastic Compton scattering (e~ + 7 — > e~ + 7) con- 
tinues to be effective until Zf w 5 x 10 4 , the photons 
can still achieve equilibrium but with a conserved pho- 
ton number. The result is a Bose-Einstein distribution 
with a non-zero chemical potential, n (rescaled by &bT 
to be dimensionless), an effect known as the "/z-type dis- 
tortion" of the black-body spectrum of the CMB, and it 
affects the distribution by 

e hv/(k B T) _ 1 ~^ e hu/(k B T)+ti _ 1 > ( 34 ) 

a positive [i reduces the number of photons at low fre- 
quencies. Finally, after Zf, even elastic Compton scat- 
tering is inefficient and photons fall out of kinetic equi- 
librium with electrons, leaving only the so-called "y-type 
distortion" [39]. As it would be difficult to distinguish 
among the y-distortions created by the heating of CMB 
photons due to diffusion damping, by the cosmic reion- 
ization (z w 10), and by the thermal Sunyaev-Zel'dovich 
effect [39] from groups and clusters of galaxies [z < 3) 
[40], we shall focus on the /i-type distortion in this paper. 

Diffusion damping occurs near the damping scale given 
as follows. Over the redshifts of interest, z ps 5x1 4 — 2 x 
10 6 , the expansion rate of the universe is dominated by 
radiation, H(z) oc (1 + z) 2 , and the effect of baryon den- 
sity on the photon-baryon fluid is negligible. Therefore, 
the damping scale, is given by [41] 

» =-/V|I±i, (35) 
J 4:5a- T n e a 45aTn e H 

which gives 

k D w 130 [(1 + z)/10 5 ] 3/2 Mpc" 1 . (36) 

Meanwhile, the heat generated by diffusion damping, Q, 
is given by 

Q = \py{fy, (37) 

where p 1 is the photon energy density and <5 7 is the pho- 
ton energy density contrast. The coefficient 1/4 mer- 
its further explanation. Naively, it would be c 2 = 1/3, 
where c s — 1/ \/3 is the sound speed of the photon fluid. 
However, a recent computation using second-order per- 
turbation theory [42] reveals that we need an additional 
factor of 3/4, yielding the number above. This heat is 
then converted into fi as 

(38) 
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The diffusion damping scales at z% and Zf are given by 
k D {zi) w 12000 Mpc -1 and k D (z f ) S3 46 Mpc" 1 , re- 
spectively. Therefore, the fi-type distortion is created by 
the (squared) photon density perturbation on very small 
scales. This property allows us to probe the power spec- 
trum on such small scales [37, 43]. 

Pajcr and Zaldarriaga recently pointed out that a non- 
zero bispectrum in the squeezed configuration makes 
the distribution of fj, on the sky anisotropic, and that 
this anisotropy of /j, is correlated with the temperature 
anisotropy of the CMB, which measurements are on 
scales much larger than the damping scale [21]. This 
allows us to measure the bispectrum in the squeezed 
configuration with a larger value of k\/k than previ- 
ously thought possible. The smallest possible wavenum- 
ber one can measure from the CMB anisotropy in the 
sky corresponds to the quadrupole, i.e., k ss 2/r^ ~ 
1.4 x 10~ 4 Mpc" 1 , where r L w 14000 Mpc is the comov- 
ing distance to the last scattering surface at zl = 1090. 
This gives k x /k S3 k D /k = 3.3 x 10 5 - 8.6 x 10 7 , which 
is far greater than that accessible from the temperature 
anisotropy of the CMB in I = 2 - 3000, i.e., h/k = 
1 — 1500, or that accessible from the scale-dependent bias 
of LSS: h/k ss l/(kR) ss 10 3 for k S3 10~ 3 h Mpc" 1 (the 
lowest wavenumber that can be plausibly measured from 
the LSS data in near future) and R S3 1 h Mpc. 



B. Cross-power spectrum of CMB temperature 
anisotropy and /ii-type distortion 



mentioned earlier, from [42]). Here, W(x) = 3ji(x)/x is 
a filter function; k s is the scale over which the damped 
acoustic waves are averaged to give heat (and which we 
will take to be equal to ko(zf) to obtain a lower bound on 
the ^-distortion); ~ 14 Gpc is the distance to the sur- 
face of last scattering; 3cos(fcr)e _fe / fe ° comes from the 
small-scale limit of the photon linear transfer function; 
and ()p denotes an average over the oscillation period. 
Correlating Eqs. (39) and (40), we find 

Cl ~20W kldkl l e \ Zf 

x J k 2 dk W (j^j B c (k u k2,k)ji{kr L )g T i(k), 

(41) 

with k s ss ko(zf). 

In order to quantify how well we can measure Cf T 
in real data, we shall estimate the cumulative signal-to- 
noise ratio, S/N, from 

(|) 2 = EWi) ^^ - (42) 

Here, we have assumed that the temperature data on 
large scales are dominated by the signal (which is already 
the case for the WMAP data), while the /u-type distortion 
data are dominated by noise. 



First, we decompose the CMB temperature anisotropy 
on the sky into spherical harmonics: 5T(n)/T = 
J2im a fm^rn{n). The spherical harmonics coefficients are 
related to the primordial curvature perturbation as 



a lm = —(-V 



d 3 k 
(2tt) e 



C(k)g Tl (k)Y; m Ck) , (39) 



where gn(k) is the radiation transfer function. Our sign 
and normalization are such that gri(k) — > — ii(fcrz,)/3 
in the Sachs- Wolfe limit. In other words, ST(n)/T — > 
— C(^ r i)/5 in the Sachs- Wolfe limit. However, we will not 
use the Sachs- Wolfe limit (except for comparison), and 
instead use gri(k) as computed from a linear Boltzmann 
code [44]. 

Next, we similarly decompose the distribution of fi 
measured on the sky, (J,(n), into spherical harmonics: 
M") = 12im a lm Y lm(n)- Following [21], we write 



^^^(fc)C(kr)C(k 2 )W (| I x 



x ji(kr L )(cos (fcir) cos (k 2 r)) p 



-{kl+ki,)/k 2 D {z) 



(40) 



with +k 2 + k = 0. (Note that the coefficient of our ex- 
pression is 187r instead of 247r because of the factor of 3/4, 



C. Estimating the noise level of /i: absolutely 
calibrated experiments 

One can relate fi to a small change in the CMB photon 
intensity, I„ = {2hv z / 'c)(e x+ » - l)" 1 , as 



2hv 3 



[1 : 



ju=0 



{e x - l) 2 



= -2.70 x 10" 18 W m- 2 Hz" 1 sr" 1 



x 3 e x 



(e*-l) 2 ' 



(43) 



where x = hv/{k B T) = 1^/(56.80 GHz), for T = 2.725 K. 
This gives 



// 



{e x - l) 2 



SI U 



10- 8 

(44) 

which can be used to estimate the noise level of /x from 
that of I v . A factor 1 x ^J is typically of order unity: 

(e *~i )2 = 1-038, 0.7307, 0.6568, and 0.8804 for v = 60, 
ldbf 150, and 240 GHz, respectively. 

For example, a proposed satellite experiment, PIXIE 
[45], is designed to have a typical noise level of 51 v — 
4 x 10~ 24 W m~ 2 Hz -1 sr _1 within each of 49152 equal- 
area pixels covering the full sky and each of 400 spectral 
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channels covering 30 GHz to 6 THz with a 15 GHz band- 
width. Averaging over the full sky, PIXIE would reach 
\fj,\ w 0.5 x 10~ 8 at 100 GHz. 

We shall assume that noise is white. The white noise 
level in the power spectrum can be calculated as: 



N 



(1(7- uncertainty in fi per pixel) 2 

x (Solid angle of a pixel in units of steradians) 

x V 2 , (45) 



where bi is the the so-called beam transfer function, 
which is the spherical harmonics coefficients of an ex- 
perimental beam profile. For a Gaussian beam with a 
full-width-at-half- maximum of 6b, h is given by 



h = exp 



Pel 

16 In 2 



(46) 



For PIXIE, we shall take the la- uncertainty in \i averaged 
over the full sky to be 10~ 8 . Therefore, the white noise 
level in the power spectrum is given by 47r x 10~ 16 [21]. 
Their beam has b = 1.6°, yielding C^' N = 4ttx 10~ 16 x 



D. Estimating the noise level of /i: relatively 
calibrated experiments 

One must have an absolutely-calibrated experiment 
such as PIXIE in order to measure a uniform /U. However, 
an interesting implication of a non-zero bispectrum in 
the squeezed configuration is that /i becomes anisotropic. 
This induces a position-dependent temperature fluctua- 
tion as hv/{kT) -> hv/{kT) + /i, i.e., 



T(h) 



2.725 K 



1 + 



The level of anisotropy is thus 
ST(h) 



S/j,(n) 



(47) 



(48) 



where Sfj, is a fluctuating part of jj,, i.e., fi(n) = p,+Sjj,(n). 
Therefore, in principle, experiments which are calibrated 
to the CMB dipole such as WMAP and Planck, as well as 
the proposed LiteBIRD [46] , are also capable of measur- 
ing this effect by making a map of S/i(n) from the differ- 
ence between temperature maps at two different frequen- 
cies, i>i and v 2 - Then, the formula for the noise power 
spectrum becomes 



N 



V\V2I{V\ - V-l) 

56.80 GHz 

x [(la- uncertainty in ST/T per pixel at vi) 2 

+ (lcr-uncertainty in ST/T per pixel at V2) 2 ] 
x (Solid angle of a pixel in units of steradians) 
x V 2 . (49) 



According to the Planck Blue Book [47] , the expected 
sensitivities of Planck at 100 and 143 GHz are 5T/T = 
2.5 x 10~ 6 per 10' x 10' pixel, and 5T/T = 2.2 x 10~ 6 per 
7.1' x 7.1' pixel, respectively. The in-flight performance 
then shows that the achieved noise level is 70% of the 
expectation, and the beam sizes at 100 and 143 GHz are 
9.4' and 7.2', respectively [48]; thus, we estimate Planck's 
sensitivity to S/j, measured from maps at 100 and 143 GHz 
as C^' N w 1.1 x 10~ 15 x e' 2 / 861 ". This is comparable 
to the above estimate for PIXIE, which is based on the 
absolute measurement of the CMB spectrum. 

However, as the sensitivity of PIXIE (and LiteBIRD) 
to CMB anisotropy is at least an order of magnitude bet- 
ter than that of Planck, if we focus only on the spatially- 
varying part of /j, rather than the uniform part of fx, then 
it may be possible to increase the sensitivity to C^ T . The 
expected noise power spectrum is of order C^' w 10~ 17 
or better, according to Eq. (49). In other words, the 
signal-to-noise of Cf can be improved by an order of 
magnitude as compared to the case where we look at ab- 
solute measurements of /i. 

In order to do this in practice, we must calibrate 
instruments at different frequencies so that they have 
the equal response to the thermal CMB. To estimate 
the required precision for calibration, let us suppose 
that the response of one instrument at v\ is different 
from that of another at v 2 by e. Then, the differ- 
ence between two maps at these frequencies will yield 
(ST/T)^) - {5T/T)(v 2 ) = e(ST/T), where ST is the 
CMB anisotropy. This residual will be confused as a sig- 

e{ST/T) which, 



56.80 GHz 



nal in Sfi, such that 8/j, 
in turn, will give a contamination of Cf 7 given by 
1(1 



) /.contamination 



2tt 



2 x 10 



-10, 



v\v 2 l(v\ - IA>) 



56.80 GHz 



1(1 + l)Cf T /2ir 
2 x 10- 10 



(50) 

As we shall show below (also see [21]), C^ T from the 
local-form bispectrum is of order 1(1 + l)|C[ lT |/(27r) w 
4 x 10~ 17 /"nl- Therefore, the required precision for the 
calibration is given by 

'vxv 2 j(v\ - v 2 )~ 



56.80 GHz 



e«2x 10- 7 /] 



NL- 



(51) 



For example, the calibration precision of the WMAP 
data is e = 2 x 10~ 3 [49], and thus the contamination 
of C^ T due to the calibration mismatch is negligible for 



/nl » 10 4 



^1^2/(^1-^2) 

56.80 GHz 



To search for a C^ T signal, 

future CMB experiments such as LiteBIRD may wish 
to place more emphasis on the relative calibration of 
their instruments. Note that the signal from a modified 
initial state will be much larger, 1(1 + l)|C[ lT |/(27r) w 

2xl °" 11 (^/0.01) (see Section VIF), so the required 
precision for the calibration can be relaxed greatly. 
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Another factor which may limit the utility of relatively 
calibrated instruments is foreground contamination. In 
principle, foreground contamination can be removed by 
using multi-frequency data, since the frequency spectrum 
of the foreground (roughly proportional to v~ 3 , v~ 2 , and 
v 2 for synchrotron, free-free, and dust emission, respec- 
tively) is different from that of fi-type distortion (oc v~ x )- 
Detailed study of foreground contamination is beyond the 
scope of this paper, but it is worth studying this in detail 
once the required calibration precision is reached. 



E. Results for the local-form bispectrum 



4 



- full transfer function 
■ Sachs-Wolfe approx. 



1110 
1 



Before we discuss Cf from the modified initial state 
effect, let us discuss Cf from the local- form bispectrum. 
Using the local- form bispectrum (Eq. (3)) in the expres- 
sion for Cf (Eq. (41)), we find 



r ^T _ 81 
25tt 6 



NL 



/ jfcjdfci 


' e -2kt/kl(z)~ 


lo 





k 2 dk w 



k D {zf) 



Pc(h) 



Pc(k)ji(kr L )g T i(k). (52) 



For a scale-invariant spectrum, Pq 
2.4 x 10~ 9 , we find 



^A 2 with A 2 



324tt 
25 



^ w 

k 



kp(zj) 
k D (zf) 
k ' 
k D (z f )_ 



ji(kr L )g T i(k). (53) 



While we use the exact radiation transfer function cal- 
culated from the linear Boltzmann code (except as indi- 
cated), it is instructive to obtain an analytical expression 
for the Sachs- Wolfe limit, gri(k) — > — ji(krL)/3. As the 
wavenumbers that are responsible for the Sachs- Wolfe 
regime are much smaller than kn(zf), one can approxi- 



mate W 



C 



k D (z 



1. We then find 



547T 

~25~ 



/nlA^Iu 



-3.5 x 10- 17 /nl x 



kp(zj) 
k D (z f ) 
2tt 



1(1 + 1) 



1(1 + 1) 



(54) 



This result agrees with that obtained by [21] up to a fac- 
tor of 3/4 recently found by [42]. Therefore, on large 
scales where the Sachs- Wolfe approximation is valid, the 
cross-power spectrum is "scale invariant," in a sense that 
1(1 + l)Cf ^constant. The overall sign is negative for a 
positive /nl because, for a positive curvature perturba- 
tion £, 5T/T = — C/5 is negative whereas the fluctuation 
in S/j, oc /nlC is positive for a positive /nl- 

How good is the Sachs- Wolfe approximation? In Fig- 
ure 3, we compare the Sachs- Wolfe approximation with 
the exact calculation. We find that the Sachs- Wolfe ap- 
proximation breaks down at I w 10, and the acoustic 



FIG. 3. Cf, the cross-power spectrum of the /i-type distor- 
tion and the CMB temperature anisotropy, from the local- 
form bispectrum with /nl = 1. The solid line shows Gf us- 
ing the full radiation transfer function, while the dot-dashed 
line shows it using the Sachs-Wolfe approximation. The am- 
plitude of Cf is linearly proportional to /nl- 



w I ft 





■ full transfer function (ideal) 
Sachs-Wolfe approx. 

■ full transfer function (Pixie) 



60 80 



100 



120 



140 



FIG. 4. Signal-to-noise ratio of Cf from the local-form bis- 
pectrum with /nl = 1- The bottom two lines show Cf using 
the full radiation transfer function, while the dot-dashed line 
shows it using the Sachs- Wolfe approximation. The solid, 
dashed line is for 9b — (ideal) and 1.6° (PIXIE), respec- 
tively. The noise level is C^' N = 4tt x i - le e i2e ?/(8in2) for 
all cases (i.e., the r.m.s. uncertainty of [i averaged over the 
full sky is 10~ 8 ). The signal-to-noise is proportional to /nl 
and is inversely proportional to the r.m.s. uncertainty of fi. 



40. As Cf 



oscillation changes the sign of Cf at / 
crosses zero at / 40, we expect the signal-to- noise ratio 
to grow more slowly with increasing multipole than with 
the Sachs- Wolfe approximation. 

In Figure 4, we show the cumulative signal-to-noise 
ratio of Cf from the local-form bispectrum as a func- 
tion of the maximum multipole, Z max . We find that the 
Sachs- Wolfe approximation overestimates the signal-to- 
noise ratio by about 40%. Due to the sign change in Cf, 
the signal-to-noise ratio does not grow between l max m 20 
and 80. 



For PIXIE's specification with 9 b = 1.6° and C, 
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4ir x 10 16 , the signal-to-noise ratio reaches S/N — 
4.3 x 10~ 4 /nl- Therefore, PIXIE would be able to see 
the signal if /nl ^ 2300. If PIXIE's detectors can be cal- 
ibrated so that the difference between maps at different 
frequencies cancels the CMB anisotropy to the required 
precision, then S/N can improve by an order of magni- 
tude. Reducing the beam size would not help much be- 
cause the signal-to-noise ratio grows only logarithmically 
with Z max [21]. 



- local form 
modified initial state 




F. Results for the modified initial state bispectrum 



We now calculate Cf from the modified initial state 
bispectrum. We start from 

27 H 6 



c; 



20tt j 



/ dki 


' e -2kl/k 2 D (z) 




/o 




z f 



00 dk 

— w 

k 2 



k 



ji{kr L )g T i{k)K(kr] ) 



27 



H 



40V2tt 5 / 2 2 
r °° dk 
k 2 



W 



k D (z f )_ 

\(N)[k D ( Zi ) - k D {z f )\ 
k 



k D {z f ) 



ji(kr L )g T i(k)K(krjo)(55) 



where £,(N) is given by the second lines of Eqs. (19) and 
(20) and K(kr]o) is given by the last lines of Eqs. (19) and 
(20). Once again, a factor H 6 /(f> 2 will be calculated from 
the normalization of P^(k) for a given e, N, and 8 (see 
Eq. (33)). 

To obtain an order-of-magnitude estimation, let us set 
k = 1 and take the Sachs- Wolfe limit. We then obtain 

H ' -Z(N)[k D (zi) - k D (z f )} 



C, 



40V2tt 5 / 2 2 
dk 

¥ 
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jf(kr L ) 



H 6 Z(N)[k D (zj) - k D (z f )]r L 

320V2vr 3 / 2 2 (I + |) (l + §) (Z-±) 



(56) 

Therefore, Cf fr om the modified initial state falls as 
Cf oc ^~ 3 , which is faster than that from the local 
form, oc I . However, the amplitude is proportional to 
ko{zi)rL ~ 1.7 x 10 8 instead oi\a[k D (zi)/k D (zf)] s» 5.5, 
which leads to a large amplification of the signal relative 
to the local form. 

In Figure 5, we compare the shapes of from the 
local-form bispectrum (solid line) and from the modified 
initial state bispectrum (dot-dashed line). As expected, 
for the low multipoles I < 40 (where C^ T is negative), 
Cf" from the modified initial state is steeper (by a factor 
of l/l) than that from the local form, whereas for high 
multipoles I > 40 (where Cf T is positive), Cf from the 
modified initial state is shallower (because it diminishes 
faster by l/l). 



FIG. 5. A comparison of the shapes of C ( mT (the cross-power 
spectra of the /i-type distortion and the CMB temperature 
anisotropy) from the local-form bispectrum (solid line) and 
the modified state bispectrum (dot-dashed line). The ampli- 
tudes of the graphs are scaled so they can both appear in the 
same figure; therefore, the overall amplitude is arbitrary. 



In order to better compare the modified initial state 
result (Eq. (56)) with the local-form result (Eq. (54)), let 
us use the power spectrum normalization, A 2 - = H 2 (l + 
27V)/(87r 2 e) (Eq. (33)), to rewrite Eq. (56) as 



9tt 5 / 2 4 £(AQ [kpjzj) - k D (z f )]r L 
Hh/2 ?£ (1 + 27V) 2 (l + §)(£+ |) (l — 5) 
-lJxlO" 11 e £(N) 2tt 



We find the ratio 



10- 2 (l + 2iV) 2 p 



CM T,„BD 4 _ Q x 1Q 5 



e(iv) 1(1+1) 



fiT. local 



(/ 



|)/nl 10- 2 (1 + 2iV) 2 I 2 



(57) 



(58) 



For single-field slow-roll inflation, /nl ~ e ~ 10 2 . 
Therefore, the ratio (Eq. (58)) is « 5 x 10 7 /L 

We can verify this ratio in a different way. Heuris- 
tically, for a modified initial state, the enhancement of 
the bispectrum in the squeezed configuration by a fac- 
tor of ki/k 3> 1 becomes a enhancement of Cf T by a 
factor of k D (zi)r L /l = k D (z l )/k C MB,i ~ 2 x 10 8 /Z, i.e., 
the ratio of the acoustic damping wavenumber to the 
wavenumber generating the multipole I, giving close to 
the ratio above. Note that the biggest contribution to 
the signal will come from small multipoles (in particular 
the quadrupolc I = 2), so the signal boost relative to the 
local form is indeed large. 

In Figure 6, we show the signal-to- noise ratio of Cf 
for a modified initial state as a function of the occupation 
number N. For PIXIE's specification with 6>b = 1.6° and 
c w,N = 4?r x 10 -i6 e ;X/(8in2) ; we find a large signa i_ to _ 

noise ratio: when the occupation number is N ~ 0.5, 
the signal-to-noise ratios are S/N ps 60(e/10 -2 ) and 
500(e/10~ 2 ) for 9 k w krj and k sa const, respectively. 
(Recall that, for the latter case, we have chosen a con- 
stant 9k so that it maximizes the signal.) This large 
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FIG. 6. Signal-to-noise ratio of Cf T from the modified initial 
state for e = 10 -2 , as a function of the occupation number 
N . The top and bottom lines are for 0k ~ const (dashed) 
and 8k ~ krjo (solid), respectively. The noise level is Cj MM ' = 
4-7T x lCP 16 e' / 8i for both cases (i.e., the r.m.s. uncertainty 
of fj, averaged over the full sky is 10 -8 , and the beam size 
is 1.6°). The signal-to-noise ratio is proportional to e and is 
inversely proportional to the r.m.s. uncertainty of /i. 

increase in the signal-to-noise ratio relative to the local- 
form bispectrum is consistent with the analytical esti- 
mate given in Eq. (58). 

Note that these calculations were done assuming ry 
was finite, i.e., we included the last lines of Eqs. (19) and 
(20). If we ignore the exponential terms in Eq. (15) by 
taking rj — > rj + ie\rj \ for large |77 |, then the signal-to- 
noise ratio goes down to about ten for 9k ~ kr]o and N 
0.5, which is still large enough for detection. Therefore, 
we should be able to detect C^ T in the PIXIE experiment 
unless N were very small. 

If the calibration of detectors at different frequencies 
meets the requirements (see Section VI D), then PIXIE 
would be able to improve its signal-to-noise ratio for de- 
tecting C^ T by an order of magnitude. Moreover, if 
Planck's calibration meets the requirement, Planck would 
be able to detect this signal. This merits further study. 

However, one caveat should be mentioned, namely that 
these results assume that the cut-off wavenumber, fe cu t, 
lies above the scales involved in the /^-distortion. If 
fccut lay within the /i-distortion scales, this model could 
produce a measurable CMB and LSS signal but have a 
smaller-than-expected C^ T . 



VII. CONCLUSION 

We have investigated phenomenological consequences 
of a modification of the initial state of quantum fluc- 
tuations generated during single-field slow-roll inflation. 
In our model, the initial state is given by a Bogoliubov 
transformation on the standard Bunch-Davies initial vac- 
uum state. A distinctive feature of this model is that the 



bispectrum of £ in the squeezed configuration - where 
one of the wavenumbers, k, is much smaller than the 
other two, i.e., fc <C fci ~ &2 _ is enhanced by a fac- 
tor of ki/k relative to the local- form bispectrum [15]. 
This enhancement generates notable effects on the scale- 
dependent bias of LSS and on the /i-type distortion of 
the black-body spectrum of CMB. 

For LSS, the scale-dependent bias goes as ki/k 3 in- 
stead of 1/k 2 , where k\ is a characteristic wavenumber 
corresponding to the short-wavelength mode in LSS for 
a given halo mass (Eq. (31)). 

For the /i-type distortion, the squeezed configuration 
bispectrum can make fj, anisotropic, which can be mea- 
sured by cross-correlating a map of /i with a map of CMB 
temperature anisotropy on large scales [21]. The mod- 
ified initial state enhances power spectrum of this 
cross-correlation by a factor of k^r^, which corrresponds 
to the ratio of the wavenumber of the acoustic damping 
scale to the wavenumber measured by CMB temperature 
anisotropy on large scales. We predict that an absolutely- 
calibrated experiment such as PIXIE can detect Cf 1 " un- 
less the occupation number is much smaller than of order 
unity. 

As this effect makes /i anisotropic, one may not even 
need an absolutely-calibrated experiment. If detectors 
at different frequencies are calibrated to have the same 
response to thermal CMB with the sufficient precision, 
then relatively-calibrated experiments such as Planck 
and LitcBIRD could detect this signal. 

We acknowledge that our derivation of the bispectrum 
from a modified initial state is limited by uncertainties 
about how to set initial conditions and how to translate 
these conditions into a proper calculational framework. 
While we think that the calculations presented in this 
paper capture plausible outcomes of a modified initial 
state, more investigation on quantum field theory with 
such a state is still necessary. That this model predicts 
such interesting signatures in LSS and the CMB moti- 
vates further study. 

Finally, while we have focused only on the bispectrum 
in the squeezed configuration in this paper, this model 
also predicts a large bispectrum in the folded limit, where 
the largest wavenumber is equal to the sum of the other 
two wavenumbers, k\ = &2 + k [26, 27, 31]. The obser- 
vational signatures that we have discussed in this paper 
should come also with the signal in the folded limit, which 
provides a powerful cross-check of the nature of the de- 
tected signal. 
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